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We present a theory of out-of-equilibrium ultrafast spin dynamics in magnetic heterostructures 
based on the s-d model of ferromagnetism. Both in the bulk and across interfaces, the exchange 
processes between the itinerant s and the localized d electrons are described by kinetic rate equations 
for electron-magnon spin-flop scattering. The principal channel for dissipation of angular momen¬ 
tum is provided by spin relaxation of the itinerant electrons. Our theory extends interfacial spin 
phenomena such as torques, pumping, and the Peltier and Seebeck effects to address laser-induced 
rapid spin dynamics, in which the effective electron temperature may approach or even exceed the 
Curie temperature. 
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Controlling spin flow in magnetic heterostructures at ultra- 
,—1 fast time scales using femtosecond (fs) laser pulses opens in- 
triguing possibilities for spintronics [1]. In these experiments 
£3 the laser-induced perturbations [2] stir up the most extreme 
regime of spin dynamics, which is governed by the highest en¬ 
ergy scale associated with magnetic order: The microscopic 

a - spin exchange that controls the ordering temperature (Tc). 
* In contrast, at microwave frequencies the ferromagnetic dy- 
namics in the bulk are well described by the Landau-Lifshitz- 

S Gilbert (LLG) phenomenology [3], which has been success- 
( fully applied to the problem of the ferromagnetic resonance 
(FMR) [4]. At finite temperatures below Tc, the spin See- 
£h beck and Peltier effects [5, 6] describe the coupled spin and 
O heat currents across interfaces in magnetic heterostructures. 
These interfacial effects are governed by thermally activated 
magnetic degrees of freedom with characteristic frequencies 
that are much higher than the FMR frequency. Despite their 
different appearances, the microwave, thermal, and ultrafast 
spin dynamics are all rooted in the exchange interactions 
m between electrons. It is thus natural to try to advance a mi- 
croscopic understanding of the ultrafast spin dynamics based 
on the established phenomena at lower energies. 

Although some attempts have been made [7, 8] to ex- 
tend the LLG phenomenology to describe ultrafast demag- 
. , netization in bulk ferromagnets, no firm connection exists 
between the ultrafast spin generation at interfaces and the 
k 4 microwave spin-transfer and spin-pumping effects [9] or the 
thermal spin Seebeck and Peltier effects. In this Letter, 
we unify the energy regimes of microwave, thermal, and ul¬ 
trafast spin dynamics in magnetic heterostructures from a 
common microscopic point of view. In the ultrafast regime, 
rapid heating of itinerant electrons leads to demagnetiza¬ 
tion of localized spins via electron-magnon spin-flop scatter¬ 
ing. The parameters that control the high and low energy 
limits of spin relaxation originate from the same electron- 
magnon interactions. In addition to the unified framework, 
this Letter’s unique contributions are the history-dependent, 
non-thermalized magnon distribution function and the cru¬ 
cial role of the out-of-equilibrium spin accumulation among 
itinerant electrons as the bottleneck that limits the dissipa¬ 
tion of spin angular momentum from the combined electronic 


system. 

The first reports on ultrafast demagnetization in Ni [10] 
challenged the conventional view of low-frequency magneti¬ 
zation dynamics at temperatures well below Tc . A multitude 
of mechanisms and scenarios have been suggested to explain 
the observed quenching of the magnetic moment. Some ad¬ 
vocate direct coherent spin transfer induced by the irradi¬ 
ating laser light as the source of demagnetization [11]. Al¬ 
ternative theories argue that ultrafast spin dynamics arise 
indirectly through incoherent heat transfer to the electron 
system [12, 13]. Recent experiments have demonstrated that 
non-local laser irradiation also induces ultrafast demagneti¬ 
zation [14], and atomistic modeling [15] supports the view 
that heating of magnetic materials is sufficient to induce ul¬ 
trafast spin dynamics. Although it is unlikely that a single 
explanation applies to all experimental scenarios of ultrafast 
demagnetization, it may be possible to address many mea¬ 
surements of different magnetic materials from a common 
microscopic view. 

THz magnon excitations in metallic ferromagnets have 
recently been proposed as an important element of ul¬ 
trafast demagnetization by several authors [16, 17]. The 
elementary interaction that describes these excitations is 
the electron-magnon scattering. Our theory is based on 
kinetic equations for the low-frequency spin and charge 
transport associated with the microwave magnetization dy¬ 
namics in heterostructures [18] and with the linear spin- 
caloritronic response [6, 19]. We extend these theories to 
treat high-temperature far-from-equilibrium spin dynamics, 
in which transport is dominated by magnons and hot elec¬ 
trons. Electron-magnon scattering plays a critical role in this 
regime. We base our understanding of this interaction on the 
transverse spin diffusion [20] in the bulk and the spin-mixing 
physics , e.g., spin transfer and spin pumping [18, 21], at the 
interfaces. 

In our approach, we assume that the localized spins that 
result in macroscopic magnetization are distinct from the 
itinerant electron bath at the energy scales of interest. In 
turn, the electron bath acts as a reservoir for absorbing and 
dissipating spin angular momentum. The experimental tech¬ 
niques that are used to measure ultrafast spin dynamics typi- 
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cally probe the absorption edge of localized electron orbitals 
buried deeply below the Fermi level [22]. Therefore, these 
methods may be less sensitive to the itinerant electron spin 
density deposited in the vicinity of the Fermi energy, such 
that the electron reservoir becomes partially invisible, or 
’’dark,” to these element-specific probes. We assume that 
the relaxation of the spin accumulation among the itinerant 
electrons is considerably slower than the spin transfer from 
the localized to the itinerant electrons. In this way we avoid 
introducing any novel ultrafast channels for dissipation of 
spin angular momentum. 

According to the accepted description of relaxation in fer¬ 
romagnetic metals, the loss of energy and angular momentum 
from localized d electrons is mediated by the exchange inter¬ 
action to the itinerant s electrons. The spin transfer from d 
to s states is accompanied by the relaxation of the s electron 
spins to the lattice through an incoherent spin-flip process 
caused by the spin-orbit coupling. Eventually, most of the 
heat stored in the itinerant electron system escapes into lat¬ 
tice vibrations (phonons). Mitchell formulated such a model 
several decades ago to describe the longitudinal relaxation 
of ferromagnetic metals [23] . A similar description was later 
employed to describe Gilbert damping in ferromagnets at low 
frequencies [24, 25]. 

We proceed by outlining the basic formalism for ferromag¬ 
netic metals in the bulk, and we will later extend the model 
to describe the spin flow across a ferromagnet (F) | normal- 
metal (N) interface, which is important for ultrafast spin dy¬ 
namics in magnetic heterostructures [1]. The Hamiltonian 
that describes F is 

H = H 0 + H sd , (1) 

where Ho consists of decoupled s- and d-electron energies, in¬ 
cluding the kinetic energy of the itinerant electron bath, the 
d-d exchange energy, dipolar interactions, and the crystalline 
and Zeeman fields. The s-d interaction is 

Hsd = Jsd E ' s ( r j) I (2) 
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where J s d is the exchange energy and [s(rj)] is the d- 
electron ( 5 -electron) spin vector (spin density) at lattice 
point j. We express the s-d interaction in terms of bosonic 
and fermionic creation and annihilation operators: 

H s d = ^ ^ Vqkk'OjqC^^Ck' 4 , T H.C. , (3) 

qkk' 

where aj (a q ) is the Holstein-Primakoff creation (annihila¬ 
tion) operator for magnons with wavenumber q and c\ a (cfc CT ) 
is the creation (annihilation) operator for s electrons with 
momentum k and spin a. H s d describes how an electron 
flips its spin while creating or annihilating a magnon with 
momentum q and spin h. The scattering strength is deter¬ 
mined by the matrix element V q kk f • 

In Eq. (3), we have disregarded terms of the form ^ 
a q a q' c ka c k'cn which describe multiple-magnon scattering and 
do not contribute to a net change in magnetization along 
the spin-quantization axis. We have also disregarded higher- 
order terms associated with the Holstein-Primakoff expan¬ 
sion. When the s-d coupling (2) is not the dominant contri¬ 
bution to H, we follow a mean-field approach and use Fermi’s 



FIG. 1. (a) Sketch of the density of s electron states in a ferro¬ 
magnetic metal with saturation spin density S. At equilibrium, 
the exchange splitting A xc shifts the bands for spin-up and spin- 
down electrons, (b) A pulsed laser heats the s electron bath. 
The out-of-equilibrium spin accumulation /jl s = S/jl^ — S/jl^ results 
from two different mechanisms: 1) electron-magnon scattering in¬ 
duces an out-of-equilibrium spin density among the s electrons, 
and 2) the mean-field exchange splitting is shifted by SA XC by the 
induced non-equilibrium spin density rid- 


Golden Rule to compute the spin transfer between the s and 
d subsystems. We assume that all relevant energy scales 
are much smaller than the Fermi energy ep = kpTp of the 
itinerant s electrons. In this limit, the electronic continuum 
remains largely degenerate, with electron-hole pairs present 
predominantly in the vicinity of the Fermi level. 

We orient the coordinate system such that the localized 
spin density points in the negative z direction at equilibrium, 
with saturation value S (in units of Ti) in the ground state. In 
the presence of a magnon density n^, the longitudinal spin 
density becomes S z = rid — S. The magnons are assumed 
to follow a quadratic dispersion relation e q = Tiu: q = eo + 
Aq 2 , where eo is the magnon gap and A parameterizes the 
stiffness of the ferromagnet. = n(e q )5 qq r defines the 

magnon distribution function n(e q ), which is related to the 
total magnon density through rid = de q V(e q )n(e q ), where 
T>(e q ) = ^Je q — eo/(47r 2 A 3 / 2 ) is the magnon density of states. 
The integral over V{e q ) is cut off at an energy corresponding 
to the bandwidth, e^ ^ ksTc, which is the magnon energy 
at the edge of the Brillouin zone. 

Because of the s-d exchange interaction (2), the itiner¬ 
ant s electrons have a finite spin density at equilibrium; see 
Fig. 1. One of the key driving forces of the out-of-equilibrium 
spin dynamics is the spin accumulation /i s = 5/jl^ — 5/ip 
The bands for spin-up and spin-down electrons are split 
by A xc ~ JsdSa 3 , where a is the lattice constant of F. 
By introducing a dynamic exchange splitting, we can write 
fji s = 5n s /D + 5A XC [26], where Sn s is the out-of-equilibrium 
spin density of the s electrons, D = 2 D^D^/(D^ + Z)^), and 
D^-m is the density of states for spin-up (spin-down) elec¬ 
trons at the Fermi level. Because the mean-field band split¬ 
ting due to the s-d exchange vanishes when the d orbitals are 
fully depolarized, 5A XC /A xc = where the sign deter¬ 

mines whether the s and d orbitals couple ferromagnetically 
(—) or antiferromagnetically (+). 
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The rate of spin-transfer (per unit volume) between the s 
and d subsystems due to electron-magnon spin-flop processes 
is determined from Eq. (3) by Fermi’s Golden Rule [21]: 

Isd = f dCqV{€q){€q /i s )X^(6g) Ms) ^l{^q)\ i 

J eo 

(4) 

where T(e q ) paremeterizes the scattering rate at energy e q . 
When the kinetic energy of the itinerant electrons equili¬ 
brates rapidly due to Coulombic scattering, ube^q — Ms) = 
{exp[/3 e ff (e q — Ms)] — 1} _1 ; this is the Bose-Einstein distribu¬ 
tion function for the spin-polarized electron-hole pairs at the 
effective temperature T e ff = 1 /(ksPett)- 

However, when the time scale of the s-d scattering is faster 
than the typical rates associated with magnon-magnon in¬ 
teractions, magnons are not internally equilibrated shortly 
after rapid heating of the electron bath, as also predicted 
by atomistic modeling [27]. Consequently, the occupation 
of the magnon states could deviate significantly from the 
thermalized Bose-Einstein distribution on the time scale of 
the demagnetization process. Our treatment of this central 
aspect differs from that of Ref. [28], in which the excited 
magnons are assumed to be instantly thermalized with an 
effective spin temperature and zero chemical potential and 
the thermally activated electron bath is assumed to be un¬ 
polarized. 

Neglecting any direct relaxation of magnons to the static 
lattice or its vibrations (i.e., phonons), d t nd = I s d/h . The 
equations of motion for the s-electron spin accumulation and 
the d-electron magnon distribution function are 

dtPs = —— —j- vCd , (5) 

T s Tl 

T f 6 j 

dtn(e q ) = - /z s ) [n B E(e 9 - /z s ) - n(e q )\ , (6) 

where p determines the feedback of the demagnetization 
on p s and r s is the spin-orbit relaxation time for the s- 
electron spin density relaxing to the lattice. r s is typi¬ 
cally on the order of picoseconds [29], and here, it repre¬ 
sents the main channel for the dissipation of angular mo¬ 
mentum out of the combined electronic system. In general, 
r s also depends on the kinetic energy of the hot electrons 
after laser-pulse excitation. This dependence will necessar¬ 
ily influence a full description of ultrafast demagnetization 
scenarios. However, this discussion is beyond the scope of 
this Letter, and we assume that r s is independent of energy. 
P — Pd -\~ Pa = — 1 /D ± A xc /S includes effects arising from 
both the out-of-equilibrium spin density and the dynamic ex¬ 
change splitting. For ferromagnetic (—) s-d coupling, these 
effects add up, whereas for antiferromagnetic (+) coupling, 
they compete. When 1 /D = A XC /S, p vanishes and p s de¬ 
couples from the dynamics of 

The s-d scattering rate can be phenomenologically ex¬ 
panded as T(e q ) = To + x{ e q ~ e o), where To (which vanishes 
in the simplest Stoner limit [20]) parameterizes the scatter¬ 
ing rate of the long-wavelength magnons and x{ e q ~~ e o) oc q 2 
describes the enhanced scattering of higher-energy magnons 
due to transverse spin diffusion [20]. In general, one might 
expect other terms of higher order in q to be present in this 
expansion as well. We will, however, limit ourselves to ex¬ 


trapolating T(e g ), which is linear in e g , up to the bandwidth 
e^, which should be sufficient for qualitative purposes. 

At low temperatures, low-frequency excitations keep the 
spin system close to a local equilibrium, resulting in purely 
transverse dynamics. In the classical picture of rigid mag¬ 
netic precession, the transverse relaxation time r 2 is de¬ 
termined by the longitudinal relaxation time t\ as follows: 
I /72 = 1/(2ti) = aw, where a is the Gilbert damping pa¬ 
rameter and uj is the precession frequency. Indeed, in the 
limit (cj,T) —* 0, Eq. (4) yields 


* _ r o 

Otn d -► — r e 0 n d , 

n 


(7) 


which is identical to the LLG phenomenology, indicating that 
eo = Tiuj and, thus, T 0 = 2a. This result establishes the im¬ 
portant link between the scattering rate T 0 in the present 
treatment and the Gilbert damping parameter that is acces¬ 
sible through FMR experiments. We note that, according to 
Eqs. (4)-(6), the itinerant electrons exert no feedback on the 
d electron dynamics in the adiabatic limit (cu, T) 0 as long 
as r s > 0. 

In the opposite high-frequency limit, pertinent to ultra¬ 
fast demagnetization experiments, we consider F to be in a 
low-temperature equilibrium state before being excited by 
a THz laser pulse at t = 0, upon which the effective tem¬ 
perature of the itinerant electron bath instantly increases 
such that T e ff > Tc . This regime is clearly beyond the 
validity of the LLG phenomenology, which is designed to 
address the low-energy extremum of magnetization dynam¬ 
ics. Dissipation in the LLG equation, including relaxation 
terms based on the stochastic Landau-Lifshitz-Bloch treat¬ 
ment [13, 30], is subject to a simple Markovian environment 
without any feedback or internal dynamics. This perspec¬ 
tive must be refined for high frequencies when no subsys¬ 
tem can be viewed as a featureless reservoir for energy and 
angular momentum. Ultrafast heating of F rapidly excites 
the itinerant electron bath to a far-from-equilibrium state. 
The subsequent buildup of p s via electron-magnon scatter¬ 
ing (4) rapidly depletes the number of available scattering 
states; see Fig. 2a. Consequently, the total relaxation of spin 
angular momentum from the combined electronic system is 
ultimately bounded by r s _1 - 

To appreciate the non-thermalized nature of the excited 
magnons, we consider the limit in which p s is small compared 
with the magnon gap eo and no magnons are excited [n(e q ) = 
0] for t < 0. After rapid heating of the itinerant electron bath 
at t = 0, the time evolution of the non-thermalized magnon 
distribution follows 


n(e q ,t) 


nBE^q, t) 


_g r( e q) e <M 


(8) 


This result implies that high-energy states are populated 
much faster than are low-energy states; see Fig. 2b for nu¬ 
merical solutions of Eqs. (5) and (6) when T e ff is increased 
from 10 2 K to 10 3 K within 50 fs and the effects of a finite 
Ms are included. By comparison, internal magnon-magnon 
interactions equilibrate the distribution function on the time 
scale ~ h~ 1 e rn [e rn /(kBTc )] 3 [21], where e m is a char¬ 
acteristic energy of the thermal magnon cloud. For short 
times, I s d (4) dominates the magnon dynamics, and we ex- 
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FIG. 2. Numerical solutions of Eqs. (5) and (6) after the itiner¬ 
ant electrons are heated from 10 2 K to 10 3 K (Tc ) within 50 fs. 
e 0 = 5 meV, A — 0.6 meV • nm 2 , p = 6 meV • nm 3 , t s = 2 ps, 
and a* = 10 a = 0.1. (a) The itinerant electron-hole pair dis¬ 
tribution n be(c — p s ) is rapidly depleted by a large spin accu¬ 
mulation p s that is built up via electron-magnon scattering, (b) 
The non-thermalized magnon distribution function n(e q ) equili¬ 
brates toward nBE(e q — p s )- However, whereas the high-energy 
magnon states are rapidly populated, the low-energy states re¬ 
main unaffected on short time scales, (c) Time evolution of the 
out-of-equilibrium spin accumulation, p s (t) : and (d) the longitu¬ 
dinal spin density, — S z (t ), after the temperature is exponentially 
relaxed toward 10 2 K on time scales of 0.15 ps, 0.5 ps, and 2 ps. 


pect the magnon population to significantly differ from the 
thermalized Bose-Einstein distribution. 

When T e ff > Tc, the thermally excited electron-hole pairs 
are populated in accordance with the classical Rayleigh- 
Jeans distribution, nsE^q ~ AG) —> ksT e ff/(e g “Ms)- As¬ 
suming, for simplicity, that the expansion for r(e g ) is valid 
throughout the Brillouin zone, Eq. (4) yields d t nd\t^o = 
Isd{ 0)/h = [r 0 + 3x(e& - e 0 )/5 ]k B T e ^S/h. Thus, the de¬ 
magnetization rate of the d orbitals is initially proportional 
to the temperature of the electron bath but is rapidly re¬ 
duced by the lack of available scattering states for high- 
energy magnons within the time scale of the demagnetization 
process. This finding conflicts with the results obtained from 
a Langevin treatment of the LLG equation [31], in which the 
magnetization relaxation rate is proportional to the temper¬ 
ature difference at all times. Figures 2c and 2d illustrate the 
time evolution of the out-of-equilibrium spin accumulation, 
/i s (£), and the longitudinal spin density, — S z (t ), for different 
relaxation times of T e ff. 

In the ultrafast regime, the electron-magnon spin-flop scat¬ 
tering is governed by the effective Gilbert damping parame¬ 
ter a* = x{ e b~ eo). Recent experimental investigations of the 
magnon relaxation rates on Co and Fe surfaces confirm that 
high-g magnons have significantly shorter lifetimes than do 
low-g magnons [17]. It is reasonable to assume that the same 
effects are also present in the bulk. The initial relaxation 
time scale in the ultrafast regime is t* ^ 

This generalizes the result of Koopmans et al. [7] for the 
ultrafast relaxation of the longitudinal magnetization to ar¬ 
bitrary a* based on the transverse spin diffusion [20]. The 
notion of magnons becomes questionable when the intrinsic 
linewidth approaches the magnon energy, which corresponds 
to a* ~ 1. Staying well below this limit, and consistent with 



FIG. 3. Sketch of a metallic ferromagnet (F) coupled to a nor¬ 
mal metal (N). In the ultrafast regime, both the rapid heating 
of s electrons in F by A T s (scenario (1)) and the heating of N 
by ATn (scenario (2)) demagnetize the d electrons in F. I s d (4) 
leads to the s electron spin accumulation p s in F, whereas U (9) 
leads to the spin accumulation p^ a t the F|N interface. Subse¬ 
quently, pn(x) diffuses into N until it vanishes due to spin-flip 
dissipation to the lattice. The additional interfacial spin current 
I s n between the s electrons in F and the itinerant electrons in N, 
due to the thermodynamic biases 5p = p s — p% and 5T = T s — Tn, 
can be described by conventional thermoelectric parameters for 
longitudinal spin-dependent transport [33]. 


Refs. [17, 20], we use a * = 0.1. For Tc = 10 3 K the initial 
relaxation time scale ~ 10 2 (Tc/T e ^) fs, which is gener¬ 
ally consistent with the demagnetization rates observed for 
ultrafast demagnetization in Fe [32]. 

Figure 3 presents a schematic illustration of an F|N in¬ 
terface. In magnetic heterostructures, and for stand-alone 
ferromagnets on a conducting substrate, the demagnetiza¬ 
tion dynamics of F are also affected by the spin accumu¬ 
lation in N, Pn(x), which can impact how non-local laser 
irradiation (e.g., the heating of N alone) induces ultrafast 
demagnetization of F [14]. By adding terms of the form 
~ Eqkk' U q kk'a q c\pck'i to H sd , where 4 t (4'i) describes the 
creation (annihilation) of an electron with spin up (down) at 
the F|N interface, the interfacial spin transfer (per unit area) 
due to electron-magnon spin-flop scattering is [21] 

re b 

h = I de q T l (e q )(e q - ii^)V(e q ) [n B E(e g - ph) ~ n i e q )] , 
d e o 

(9) 

where p^ = /xn(0 ) is the spin accumulation at the interface 
and T l (e q ) parameterizes the interfacial scattering rate. 

Importantly, we note that the scattering of coherent long- 
wavelength magnons at the F|N interface can be described in 
the language of spin pumping/spin Seebeck effects [21] and 
can be parameterized in terms of the spin-mixing conduc¬ 
tance (per unit area) [18]. Motivated by T(e g ) in the 
bulk, we write r 2 (e g ) = g^(e q )/(nS), where g^ reduces to 
g^l for low-energy scattering, e q 0. The interface scat¬ 
tering (9) dominates the microwave spin relaxation in thin 
ferromagnetic layers of thickness dp ^ 10 nm [18, 34]. This 
trend should continue for higher frequencies and is relevant 
for ultrafast spin dynamics in thin magnetic layers in het¬ 
erostructures [1]. We expect the energy dependence of the 
effective interfacial spin-mixing conductance to be relatively 
weak compared to that of the bulk scattering T(e q ), which 
















5 


can be severely constrained at low energies due to momentum 
conservation [20]. For a finite temperature bias ST across 
the interface and for magnons thermalized at the tempera¬ 
ture T < Tcj the connection to the thermal spin Seebeck 
and Peltier effects is made by identifying <S = drh and 
II = TS/h [19] as the Seebeck and Peltier coefficients, re¬ 
spectively. In the ultrafast regime, the relative importance 
of the bulk scattering, parameterized by a *, and the interfa¬ 
cial scattering, parameterized by g^, can be extracted from 
measurements of demagnetization strength and spin currents 
in magnetic heterostructures. 

In conclusion, we have extended the concepts of trans¬ 
verse spin diffusion in bulk ferromagnets and the spin-mixing 
physics across interfaces to address the ultrafast spin dynam¬ 
ics observed in rapidly heated magnetic heterostructures. 
In the microwave, thermal, and ultrafast regimes the ele¬ 
mentary interactions that defines the demagnetization time 
scale are described by kinetic rate equations for electron- 
magnon spin-flop scattering. For metallic ferromagnets in 
the bulk, our analysis shows that treating the itinerant and 
localized electron subsystems as being individually thermal¬ 
ized with equilibrium distribution functions parameterized 
by an effective temperature is insufficient to describe the 
far-from-equilibrium spin dynamics that arise from pulsed 
laser heating. The magnon distribution function remains 
non-thermalized on the relevant time scale of the demag¬ 
netization process. We emphasize the critical role of the 
out-of-equilibrium spin accumulation g s among the itiner¬ 
ant electrons, which provides the bottleneck that limits the 
relaxation of spin angular momentum from the combined 
electronic system. 

The authors thank A. V. Kimel, G. E. W. Bauer, J. Barker, 

S. Bender, and H. Skarsvag for valuable discussions. 
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